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We previously established that in principle, it is possible to quantum compute using passive linear optics with 
photo-detectors [JlJ. Here we describe techniques based on error detection and correction that greatly improve 
the resource and device reliability requirements needed for scalability. The resource requirements are analyzed 
for ideal linear optics quantum computation (LOQC). The coding methods can be integrated both with loss de- 
tection and phase error-correction to deal with the primary relaxation processes in non-ideal optics, including 
detector inefficiencies. The main conclusion of our work is that the resource requirements for implementing 
quantum communication or computation with LOQC are reasonable. Furthermore, this work clearly demon- 
strates how special knowledge of the error behavior can be exploited for greatly improving the fault tolerance 
and overheads of a physical quantum computer. 



I. INTRODUCTION 

In we proposed that linear optics quantum computation 
(LOQC) is a viable option for physically realizing quantum 
computers. The proposal depends on a series of optical pro- 
tocols that require single photon state preparation and mea- 
surement whose outcome can be used to control other opti- 
cal elements. The basic idea is to kick back the hidden non- 
linearities in photo-detectors to qubits encoded in pairs of op- 
tical modes (photonic qubits). The protocols were shown to 
implement the necessary quantum gates with arbitrarily high 
probability of success. In their simplest form, the resources 
required for implementing the protocols grow rapidly as the 
desired success probability is increased. We noted that due to 
the general accuracy threshold theorem [^, ||, g], the model 
scales efficiently asymptotically, with constant overhead for 
implementing the standard fault tolerant model of quantum 
computation. We also suggested that with the use of erasure 
codes [||, this overhead could be significantly reduced. 

Here we apply the techniques of of quantum error-detection 
and correction to greatly increase the efficieny of implement- 
ing quantum information processing by LOQC. We first show 
that for ideal LOQC (where errors in devices are ignored), it 
is possible to use a two qubit error-detecting code to rapidly 
eliminate the probability of failed implementations of the two 
qubit gates. This follows from a threshold analysis demon- 
strating a threshold of .5 for an error model where the errors 
are a z measurements at known locations. One consequence of 
this result is that the non-deterministic gates of [|l]| need only 
be implemented using the 2- or 3-photon prepared states. The 
next task is to demonstrate that the main sources of device 
errors can be efficiently eliminated. We accomplish this by 
adding phase error-correction, and more importantly erasure 
coding methods. The latter serves the purpose of removing 
errors caused by particle loss and detector inefficiency, ex- 
ploiting the built-in loss detection capabilities of basic LOQC. 
In the process we give a conservative bound on the threshold 



* knill@lanl.gov . http : / /www . c3 . lanl . gov : ~ knill 



I .afiamme@lanl.i 



nilbum@physics.uq.edu. an 



for the erasure error model. 

On the basis of this paper and [[jj we can propose a roadmap 
for experimental and theoretical work toward implementing 
LOQC. It is based on viewing basic LOQC as a fundamental 
model of computation that can be used to implement standard 
quantum computation by using a set of layered techniques. 
For benchmarking purposes, the relevant resources of LOQC 
can be taken to be the number of particles independently gen- 
erated, the probability of (detected) failure of the implemented 
computation (without post-selection), and the measured error 
in the output conditional on success (i.e. with post-selection). 
The roadmap for LOQC based quantum computation can be 
outlined as follows: 

1. Basic LOQC. 

a. Non-deterministic non-linear sign changes. 

b. Preparation of the state \t n ) for teleportation 
(see [[!]] for the definition of \t n )). 

c. Teleportation using \t n ) for increasing n with 
probability of success close to l/(n + 1). 

d. Controlled sign changes with good probability of 
success. 

2. QC based on LOQC. 

a. Boosting the success probability by the use of en- 
coding. 

b. Decrease phase errors by applying error correc- 
tion. 

c. Decrease loss errors by encoding with erasure 
codes. 

d. Concatenation (or larger codes) for achieving high 
accuracy. 

The experimental challenges include the ability to use mul- 
tiple independently generated single photons and to control 
the photonic qubits using feedback from detectors. The first 
demonstrations are likely not to involve feedback, but rather 
to use post-selection and high repetition rates to demonstrate 
success. Although feedback can be delayed in our schemes, 
this comes at the cost of high failure rates in the state prepa- 
ration protocols, particularly for those states that depend on 



2 



previously prepared states. The ultimate goal is to build the 
state preparation protocols into state factories with the ability 
to attempt failure-prone state preparations at high rates (per- 
haps in parallel), exploiting the ease with which photons can 
be generated. 

Here is the outline of this paper. We begin by recalling 
the needed properties of the LOQC protocols given in [[!]]. 
We briefly describe the concatenation method for establish- 
ing thresholds, give a code suitable for increasing the success 
probability in ideal LOQC, and show how to implement en- 
coded operations. The gain in success probability is estimated 
so as to determine a threshold. The resources that determine 
the general error propagation behavior are bounded for the 
goal of approaching the quantum communication threshold. 
As a result we conservatively estimate that a gate with the 
necessary reliability depends on at most a few hundred LOQC 
controlled sign flips implemented non-determinstically with a 
pair of states each involving three photons. More are used up 
in unsuccessful state preparations, but at least in principle, this 
overhead is not much larger. The quantum coding methods 
are then enhanced for the purpose of dealing with phase and 
detected-loss errors, including particle detector inefficiencies, 
and finally for dealing with any residual general errors. 

We assume familiarity with quantum computation [^|, |[| 
and quantum error-correction via stabilizer codes [^, For 
the basic ideas of LOQC see [[!]]. Most of this paper is written 
in the language of qubits using products of Pauli operators. 

n. FEATURES OF LOQC 

LOQC is a model of quantum computation where some of 
the gates are non-deterministic, with detectable failures. The 
probability of failure of the gates depends on the resources 
used. Specifically, the model is characterized by enabling the 
following operations on standard qubits (which are encoded 
in the physical system as photonic qubits): 

1. Preparation of 1 0) : Po- 

2. Measurement in the basis |0), |1): Mg. 

3. Every one qubit rotation. 

4. Controlled sign flip: c-a z , with a probability of failure, 
which is detected. 

In practice errors other than detected failures occur. For the 
moment, we assume that such errors are significantly less 
likely than detected failures, so that initially, they can be ig- 
nored. Specifically, this leads to designing implementations 
by dealing with errors in order of their importance. We call 
LOQC with only the errors due to detected failure of the con- 
ditional sign flip "ideal LOQC" (iLOQC). In iLOQC, single 
photon state preparation, particle number detection and pas- 
sive linear optical elements are all perfect. As explained in 
our previous paper when applying c-aj 12 \ with prob- 
ability /, qubit 1 is measured in the a z basis. If this event 
does not occur then with probability /, qubit 2 is measured in 
the a z basis. Thus the prior probability of the second event 



is /(l — /). The measurement outcome is known in either 
case and the qubit not measured is preserved. We take that to 
be the error model of iLOQC. Note that the failure behavior 
is asymmetric, so that the ordering of the labels is significant. 
We call the first qubit in this operation the "source", and the 
second the "target". The resource overhead for implementing 
this gate depends on /. If the methods of [[!]] are used for 
implementing the c-a z operation, a state with approximately 
(1//) — 1 photons needs to be adjoined. The preparation of 
the state needs to be tried several times, using ancillary pho- 
tons. With the naive method, the total number of photons used 
in the preparation attempts grows exponentially in 1//, so it is 
desirable to show that we can scale with as high a probability 
of failure as possible. 

For the present purposes, it is convenient to use a gate set 
based on the product operator formalism [[Tl|], Thus, gates are 
exponentials of products of Pauli operators. To simplify the 
notation, define X := a x ,Y := (T y ,Z := a z . For a prod- 
uct operator U, write Ue = e~ ,c/7r9 / 360 and note that since 
U 2 =I,U e = cos(tt6>/360) - i sin(7r6>/360)J7. For example, 
Xi8o is a bit flip up to a global phase. When there is no pos- 
sibility for confusion, we abbreviate UV = U^V^'K where 
the parenthesized superscripts are system labels. 

One reason for why product operators are convenient is 
because it is straightforward to follow their evolution under 
90° rotations by using appropriate triples of axes. For ex- 
ample, a Ygo rotation takes the Z axis to the X axis, so that 
(Y 90 )Z(Y- 90 ) = X. This implies that if in a quantum net- 
work, a Z measurement or a Z rotation precedes a Ygo rota- 
tion, then this is equivalent to an X measurement or rotation 
(by the same angle) after the Ygo. A complete set of gates that 
generates the determinant 1 unitary matrices is given by 

1. X rotations: X^ for any angle <\>. 

2. Z rotations: Zgg. 

3. ZZ rotations: {Z^Z^) go . 

(With the use of ancillas, Zgo may be eliminated from the set 
without loss of completeness.) The third gate is related to c-a z 
by 

(Z^Z^) 9Q = e^Z 90 ^Z 9a ^c-aJ 12 l (1) 

As a result, the ZZ rotation is readily implemented in LOQC 
up to global phases. Furthermore, we can modify the im- 
plementation to achieve the following failure behavior for 
(zMz (2) ) 9 o: With probability /, qubit 1 is measured in Z 
and qubit 2 is untouched. With probability (1 — /)/, qubit 
2 is measured in Z. If the outcome is 0, qubit 1 experienced 
a Zgg and if it is 1, a Z-gg. Note that due to the availabil- 
ity of Xgg rotations, we can also use (U^V^)±gg rotations 
for U and V either Z or Y with similar error behavior, where 
the measurements commute with the rotation. Similarly, it 
is straightforward to implement the Y±go rotations, Y mea- 
surement and Y eigenstate preparation. We will find that the 
encoded Z rotation also has a probability of failure, where a 
Z measurement occurs if it fails. In this case, X eigenstate 
preparation (using |0) preparation and an Xgg followed by a 
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Z±go) may fail. However, since the preparation is successful 
if failure has not been detected, it is possible to retry it until 
success is achieved. This makes it possible to get a perfect X 
eigenstate preparation. 

m. ESTABLISHING ACCURACY THRESHOLDS 

For the very special error model introduced in the previous 
section, accuracy threshold analyses are much simpler. The 
basic principle is to use a quantum code that permits imple- 
menting the basic operations on the encoded information in 
such a way that the new error model is consistent with the 
original one, and so that the new error rate is substantially 
less. The encoded qubits then behave just like the more fun- 
damental ones used to encode them, so that the same coding 
method can be used with these new types of qubits. This re- 
cursive coding method is known as concatenation and has the 
property that the error-rates decrease super-exponentially with 
the number of levels of concatenation. The basic components 
of an accuracy threshold result by concatenation are the fol- 
lowing: 

1 . A quantum code. 

2. Means of implementing each of the basic operations on 
the encoded qubit. 

3. Means of recovering from error, which may be part of 
2. 

4. Establishing the error model that applies to the encoded 
qubits and calculating a bound on the new error rate. 

The plan is to show that a two qubit code suffices for es- 
tablishing a threshold of / = .5 for iLOQC. After dealing 
with the errors of iLOQC, it is necessary to consider the con- 
tributions of other sources of noise, particularly photon loss 
and phase error. It turns out that the same family of codes 
can be used with phase error correction. We recall that LOQC 
comes with an effective leakage detection scheme and observe 
that the basic teleportation protocol used in LOQC can be en- 
hanced to allow detection of loss from particle detector ineffi- 
ciency. As a result, it possible to use erasure code to eliminate 
errors from photon loss. We point out that due to the special 
nature of the erasure error model, good accuracy thresholds 
apply and error rates can be readily improved with relatively 
simple codes and few levels of encoding. 

IV. A TWO QUBIT QUANTUM CODE 

Let|+) := (|0) + |l))/V2and|-) := (|0) - |1))/V2 be the 
eigenstates of X. Up to overall scale factors, the associated 
projection operators are / ± X. We continue to omit these 
scale factors in identities. To improve the failure probabilities, 
we use a two qubit quantum code with the encoding 

10} - |0x> 

= (!++> + 1— 



= (|00) + |11»/V2 (2) 

= »/V2 

= (|01> + |10»/V2, (3) 

and show how to implement the necessary operations on the 
encoded states. The encoded states define the state space of 
the "logical qubit". In the language of stabilizer codes, the 
code is defined by the stabilizer X^X^, and accordingly, 
the projection onto the code space is given by / + X^X^ 2 \ 
We use L as the label for logical (encoded) operators and 
states. With this encoding, logical operators are given by 

I< 1 » = X& = L X® (4) 

Z (D = Z (1) Z (2) =l _y(l)y(l) (5) 

yW = yCD^C 2 ) = L zWyW, (6) 

where we introduced the notation =l to denote identity when 
restricted to the code. 

For the purposes of establishing a threshold, we use the fol- 
lowing set of basic operations and assumptions: 

1. X, Y and Z eigenstate (eigenvalue 1) preparation: 

Px,Py,Pz- 

2. Y and Z measurements: My and Mz- Mz{s) means 
that the result of the Z measurement was eigenvalue s. 

3. Xiao, Yiso and Ziso rotations. 

4. Xfj, rotations. 

5. Zgo rotation, with failure probability / after the first 
encoding. 

6. (Z^ 1 ' Z^ 2 >)go rotation, with failure probability / in our 
error model. 

Operations 1. to 4. are error-free in iLOQC. The Zgo rotation 
is error free in iLOQC, but as we will see, it may fail with 
probability / with a Z measurement after the first level of 
encoding. Note that Ygo, (YZ)gg and (yy)go rotations can 
be implemented by conjugation with error-free X rotations. 

A. State preparation 

To encode an arbitrary state \ip\ to \ip} L , one can use the 
sequence of operations given by 

= YWg (Z^Y^)g a Y^-g P Z W. (7) 

To see that this works, follow the effect of the unitary gates 
on the initial operators X^\ Z^ 1 ' (basic operators associated 
with the state to be encoded) and I + Z^ (the projection onto 
the prepared state). It can be seen that A'W -> JSfW = L 
z« -» -y«y( 2 ) = L ZW and I + Z& -» I + 
X^X^ (the projection onto the code). 
The sequence 

£7(12) 

can be used to prepare encoded eigen- 
states of X, Y or Z by preparing the eigenstate on qubit 1 first. 
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The process fails with probability / + (1 — /)/ = f(2 — /), 
and can be repeated an expected 1/(1 — /(2 — /)) times to 
successfully prepare it. Later, it will be the case that the Y 
rotation in the sequence can fail, which changes the failure 
probability to /(l + (1 - /) + (1 - /) 2 + (1 - /) 3 + (1 - J) 4 ). 
The resource usage can be improved by reusing qubits not af- 
fected by the measurement in the failure. 



B. Measurement 

To measure the logical qubit in the logical basis, measure 
both of the supporting qubits. A Z measurement on both 
yields a logical Z measurement via the total parity of the two 
outcomes. Similarly, a Y measurement on the first qubit and a 
Z measurement on the second gives a logical Y measurement. 
A logical X measurement is accomplished by measuring X 
on the first qubit. All of these measurements are without er- 
ror. 



C. Logical qubit rotations 

The logical 180° rotations are implemented by applying ba- 
sic ones to each qubit. Thus 



A 180 — A 1 '180 

z (1 W 

180 — * 180^ 180- 



6 180 — £ 180^ 180 



(8) 
(9) 
(10) 



These are error-free. The ability of implementing 180's in 
this way is generic for stabilizer codes [|To|]. The rotations 
X^'^ are obtained by applying X^ $ and are error-free. To 
implement Z^ 99 , apply (Z^ Z^) 9 q. This is not error-free, 
and we show later how to use recovery from Z-measurement 
to get a smaller probability of failure. 



D. ZZ rotation 

If the logical qubits are encoded in qubits 1,2 and 3,4, re- 
spectively, the logical ZZ 90° rotation can be obtained by 
applying (Z^ Z^) 90 . This can be done by the se- 
quence 



(Z (il) 2 (£a) )oo = {Y^Z^) 9Q {Y^Z^)_ 9Q 

(zWzW) 90 

{Y^Z^) 9Q {Y^Z^). 9Q , 



(11) 



where we use the convention that the order of application is 
right to left within a line and top to bottom for multiple lines. 
Unfortunately, this does not readily yield a logical gate with 
significantly less error. To do that requires using the telepor- 
tation techniques of [10 12] ]. 



V. ROBUST TELEPORTATION 

A. Basic teleportation 

The basic teleportation protocol transfers an arbitrary state 
from qubit 1 to qubit 3 by first preparing a state on qubits 2 and 
3, then making a measurement on qubits 1 and 2, and finally 
correcting qubit 3 by applying one of the 180° rotations. Here 
is a sequence T^ 123 ) for a variant of the usual protocol that has 
better error behavior for our purposes. 



(12) 



T (123) = p z (2)p y (3) 

{Y^Z^U 

(zWyW) 90 

M z W( Si )My (1) (s 2 ) 
(U(s 1 , S2 f ) ) 180 

The source in the second coupling evolution is chosen to be 
qubit 2. The necessary correction U(si , s 2 ) can be derived by 
determining the effect of the process on an input operator. Us- 
ing the projection operators I ±Z and I ±Y for the prepared 
states and for the effects of the measurement the transforma- 
tion for an initial Z^ operator is 

Z (i)( / + Z (2) )(/ + y(3) ) ^ (I + S2 Y^)(I + Sl Z^) 
xZ^(I - Z^Z^Z^)(I - y(2) x (3)) 
x(I + S1 ZW)(I + s 2 YM). (13) 

Using the rules 

(/ + Z (a) U {b) )(I + sZ (a) ) = (I + sZ {a) ){I + sU {b) ) (14) 
(/ + sZ {a) ){I + Y (a) U {b) ){I + sZ {a) ) = (I + sZ (a) ) (15) 

and their variations (continuing to omit constants in the iden- 
tities), this evaluates to -siZ^)(I + s 2 Y^)(I + siZ^). 
Similarly, 

x {1 \i + z {2) ){i + y (3) ) (j + s 2 y {1) ){i + Sl z {2) ) 

x y(i)y(2)(j _ Z ^Z^Z^)(I - yWxO)) 
x(I + Sl Z^)(I + s 2 Y^), (16) 

which evaluates to to -s 2 X^ {I + s 2 Y^){I + Sl Z^). This 
implies that 



£7(1,1) 
£7(1,-1) 
C/(-l,l) 
U(-l,-l) 



= Y 

= X 

= Z 

= I 



(17) 



For a nice group theoretic treatment of teleportation, see Jl3|]. 

The state obtained on qubits 2 and 3 before the last rota- 
tion and measurements is a prepared entanglement denoted 
by \t e \ 3 . The idea is to use the built in error-detection and 
multiple tries to obtain the state without error. 

In order to analyze the propagation of errors, we need to 
understand the effect of an unintended Z^ or Y^ measure- 
ment before the protocol's end. Both of these commute with 
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the two rotations in the protocol. A Y ^ measurement implies 
that the net effect of the applied rotations on the third qubit 
is a Z^±qq (the sign depends on the measurement outcome). 
Nothing happens to the first qubit. If a Z^ measurement hap- 
pens, this directly applies to the first qubit. The second qubit 
experiences a Y^ 2 '±go rotation. To simplify matters we inten- 
tionally perform the Y measurement on this qubit to return to 
the first case as far as qubit 3 is concerned. 



B. Logical ZZ rotations by teleportation 

One method for implementing the logical ZZ rotation 
on qubits encoded in qubits 1,2 and 3,4 respectively is to 
first teleport the four qubits, then apply the rotation Rz = 
(ZWzWzWzW)go. Since the final step of the teleporta- 
tion protocol involves a number of 180° rotations, and Rz is 
in the normalizer of the Pauli group, we can instead apply Rz 
to the destination qubits in the four copies of \t e ) and apply 
appropriately modified 180° rotations after the teleportation 
measurement. Actually, it is better to apply Rz first, use the 
original 180° corrections and note that the overall effect is 
equivalent to Rz or R z after teleportation. Which one actu- 
ally occurred can be determined from the measurement out- 
comes. To go from R z to Rz it is sufficient to apply Z lso 's 
to each qubit. Because of the ability to retry the state prepa- 
ration until it succeeds, this reduces the problem of reliably 
implementing Rz to the problem of reliable teleportation. 



C. Error recovery by teleportation 

Our methods are designed so that the only error from which 
it is necessary to recover is a Z measurement with known out- 
come on one of the qubits. It is desirable to implement the re- 
covery so that at worst, a Z measurement occurs on the logical 
qubit. By symmetry, it is sufficient to consider a Z measure- 
ment on qubit 1. The effect of the measurement is to apply 
(/ + sZ^- 1 '), where the sign s is known. One way of restoring 
the encoded qubit is to notice that 

(I+sZW)(I+xWxW) = (I+siYWxW){I+xWxW). 

(18) 

Since the first operator is a 90° rotation around Y^X^ 2 \ it 
can be undone by applying its inverse. 

A method more easily made reliable is based on the syn- 
drome measurement technique of error correction. From this 
perspective, the unintended Z measurement creates a super- 
position of states with two syndromes, that is eigenvalues of 
X^X^l The encoded qubit can be recovered by measur- 
ing S^ 12 ) = X^X^ and if the eigenvalue is —1, applying 
Z^> . To measure S^ 12 ) we again use teleportation, reducing 
the measurement problem to a state preparation and teleporta- 
tion problem (see [[mJ] for similar ideas used to solve the more 
difficult problem of correcting unknown errors). The idea is to 
measure S^ 12 ) on the destination qubits of two copies of \t e ) 
before completing the protocol, and then infer the eigenvalue 



from the combination of all measurement outcomes. The cor- 
rection operations U of the protocol are unchanged. To see 
how this works, implement the protocol by teleporting qubit 1 
with \t e ), i6 and qubit 2 with \t e \ 7 , so that the teleported state 
ends up in qubits 6 and 7. The measurement on qubits 6 and 
7 can be implemented by the sequence 

(Y^X^)^ Mz {6 Hs)(Y^X^) 90 . (19) 

Note that the two X operators on qubit 7 occurring in the 
two-qubit rotations can be obtained by conjugating Z oper- 
ators by a Ygo rotation. If the measurement of XX results 
in s = — 1, we correct the state by applying Y^i^qZ^ 3 ' iso- 
(To see that this correction works, examine the quantum net- 
work, moving the correction operator back to the beginning 
by appropriately changing orientations of rotations by anti- 
commuting operators, and then absorbing it at the commuting 
state preparation steps.) Again, we can retry this state prepa- 
ration until it succeeds. The prepared state is now given by 
(/ + i S , ( 67 ) ) |t e ^3 6 \t e ) i7 , which is already a logical qubit state 
on bits 6 and 7. Suppose teleportation concludes successfully, 
with correction operators U(s±, S2)^ 6 ' > and C/(s3, S4)' 7 '. The 
resulting state is the same as if we had applied (I + s'S^) 
after the protocol, where s' depends on the Si (180° rotations 
only change the sign of Pauli operators when changing the or- 
der of events in a sequence). If the state to be teleported is 
in the code, then s' = —1 is impossible (as the projection is 
otherwise orthogonal to the state, resulting in a zero probabil- 
ity event), while if it is in the — 1 eigenvalue space of s^, 
then s' = 1 is impossible. Because of the relationship of s' 
to the Si, which event occurred can be determined from the 
combination of s^'s that resulted from the teleportation mea- 
surement, and the necessary correction can be applied. 

VI. ERROR ANALYSIS 

We begin by considering errors that occur in the encoded 
operations and how one should respond to such errors. 

A. Errors in recovery 

The recovery procedure is applied when an unintended Z 
measurement occurs. Suppose this occurs at qubit 1. Since 
the state of qubit 1 is now known, we can take advantage of 
this to prepare a state with the first teleportation step (that is 
the one involving qubit 1) already completed. This can be 
done error- free, given a number of attempts. Specifically, the 
first teleportation protocol can be replaced by preparing the 
destination qubit in the appropriate eigenstate of Z, and then 
applying the XX measurement to this and the target of the 
second teleportation before completing the latter. As before, 
one can use either outcome of the XX measurement, in this 
case by applying only a Zigo correction, if necessary. 

If the teleportation of the second qubit fails at the source of 
the relevant rotation, the second qubit is untouched, and we try 
again. If it fails at the target, then the second qubit is measured 
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in Z, which implies that the logical qubit is measured in Z. 
Let F r be the probability of failure. By following the different 
possible outcomes in the attempts, we get 

F r = fF r + (l-f)f (20) 

.'. F r = /. (21) 



D. The threshold 

The threshold Td for obtaining an improvement in the fail- 
ure parameter can be determined by solving Fzz = f, which 
gives 

T d = .5. (26) 



B. Errors in the implementation of c 



VII. RESOURCE ANALYSIS 



When applying (Z^ Z^)go to the qubits, the following 
can happen: 1. The first qubit is measured in the Z basis. 
In this case, apply the recovery procedure and if it succeeds, 
attempt the operation again. If it fails, the logical qubit is mea- 
sured in Z. 2. The second qubit is measured in the Z basis 
with outcome s and the first qubit experiences a Z^ s go- The 
effect is the same as if a Z^ l \q had been applied before the 
measurement, so if the subsequent recovery procedure suc- 
ceeds, then the desired operation has been applied. Otherwise, 
the logical qubit has been measured in Z. The probability of 
failing once case 2 is entered is /. The probability of entering 
case 1 is /. By following the re-attempts going through case 
1, we obtain the equation for the probability of failure 



F z 
F z 



(1- 
f(2 



f)f + f 
-/)/(!- 



+ f(l-f)F z 
/(I-/))- 



(22) 
(23) 



(£2)1 



C. Errors in the implementation of (Z^ 1 ' Z 



To avoid having to retry the operation we modify the proto- 
col for implementing the logical ZZ rotation slightly. Instead 
of preparing 4 copies of \t e ) with the ZZZZ rotation already 
applied, prepare k > 4 such copies. In the end, the destination 
qubits of the unused copies are measured in Z, so that the ef- 
fective applied rotation becomes {Z^ Ll ^ Z( L ^)±go, where the 
sign depends on the measurement outcomes. Compensate for 
a minus sign by applying Ziso's to each qubit. Let the qubits 
be encoded in qubits 1, 2 and 3, 4 respectively. Perform the 
teleportation protocols for the qubits in this order. When a 
protocol fails at the source of the critical rotation, try again 
with the next available pair of qubits in the prepared state. 
When the protocol fails at the target, the procedure depends 
on whether it is the first or second member of a pair of encod- 
ing qubits. If it is the first, attempt recovery using the second 
qubit as usual. If it is the second, do the same, but using the 
already teleported qubit as the first member. In this case, we 
do not need to re-attempt teleportation for implementing the 
rotation, as in the case of the logical Z rotation. Computing 
the probability of failure that the procedure fails by the first 
logical qubit being measured gives 

Fzz = f 2 + (l-f)f 2 + f(l-f)F ZZ (24) 
.-. F zz = F z . (25) 

The probabilities for the second logical qubit being measured 
given successful completion of the first two steps is the same, 
as required by the assumptions of the model. 



Resource analysis can be used to estimate the effect of 
residual errors (not fitting the error model) in the basic op- 
erations and to determine the total overhead of implementing 
an accurate quantum gate for standard quantum computation 
or communication. The total overhead includes the expected 
number of attempts required to prepare the requisite states. 
In an actual system, the state preparation attempts can be ar- 
bitrarily parallelized and implemented in independent high- 
throughput state factories. In the system as proposed here, 
the states are relatively simple in terms of the lower level im- 
plementation, and success probabilities are reasonable. An 
explicit total resource analysis is left as a problem for future 
work. For now, our primary concern is how general errors 
can propagate from the physical implementation to the en- 
coded qubits. This depends only on the operations that di- 
rectly contribute toward the state used in the final gate via their 
errors conditional on success. This property can be exploited 
to la rgely el iminate the problem of inefficient detectors, see 



Sect. VIII B 



The analysis that follows is intended as an example for how 
this can be done and is completed with an explicit example. 
We begin by counting resources in terms of the operations of 
iLOQC, counting separately first the error-free one qubit rota- 
tions, state preparations and measurements in the Z or Y basis 
CR ), second the 90° Z and Y rotations (R x ), and third the 90° 
two qubit rotations with Y or Z operators (i?2)- This separa- 
tion helps with the resource estimate due to the fact that they 
differ in resource requirements at the first and later levels. For 
our purposes two or three levels are expected to suffice. Note 
that we are not counting steps that are required to temporar- 
ily store a qubit. This is necessary in principle, as imperfect 
memories without parallelism imply that truly scalable quan- 
tum or classical computing is impossible [ 15 1. In particular, it 



is beneficial to parallelize implementations as much as possi- 
ble. 

As we proceed, we will comment on the expected number 
of tries for state preparations. For this purpose, define q as the 
probability of failure of the one qubit rotations contributing to 
Ri and let p = 1 — f(2 — f) be the total failure probability of 
the two qubit rotations contributing to i?2- In a total resource 
analysis, one can exploit the fact that q = at the first level. 



A. Resources for teleportation 



The preparation of \t e ) requires 

Ro{t e ) = 2 



(27) 
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Rl(te) 
R 2 (t e ) 



(28) 

(29) 



Since the probability of success is (1 — p), the expected num- 
ber of attempts to assure success is 1/(1 — p). 

To prepare 2k copies of \t e ) with the (ZZZZ . . .) 90 rota- 



tion applied to the targets using the method given in Sect. [VP 
requires 



Rv(z 2k ) 

R 1 (Z 2k ) 
R 2 (Z 2k ) 



2kR (t e 
2kR 1 {t e 
2kR 2 (t e 



= 4fc 
= 
+ 4fc- 



1 = 6fc + 1. 



(30) 
(31) 
(32) 



The probability of successful preparation is only (1 — p) 4fc+1 , 
so the preparation method needs to be improved. An efficient 
(in fc) scheme is based on the idea of using a parity containing 
ancilla to kick back the desired rotation, and to generate this 
ancilla in a tree like fashion |16|. The sequence is defined 
recursively by: S\^ ab ' consists of preparing \t e ) a and |0^, then 
applying c-a x from the target of \t e ) a to |0^. The qubit b is 
the "parity" qubit. 5 i+ i (abd) applies c-cr T (M) to the outputs of 
S , / ab - ) and S^ cd \ then measures parity qubit b in the X basis 
by applying Y 9 q and then measuring Z. If the outcome is — 1, 
apply a -Ziso to each of the target qubits of the \t e ) that make 
up a. The new parity qubit is d. The controlled-not operation 
c-<j x ( ab } is applied with 



Z^_ 90 X^ 9 oY^- 90 (Z^Z^) 9Q Y^ 



(33) 



To kick back the desired rotation to 2 l copies of \t e ) after Si 
has been successfully completed, apply Z 9 q to the parity qubit 
and measure it in the X basis, applying a Ziso correction as 
before, if necessary. The last step can be deferred until after 
the teleportation when using this for implementing the logical 
ZZ operation. The failure response of the algorithm can be 
optimized by recovering states as much as possible. For sim- 
plicity, we assume that the state associated with a failed Si is 
discarded. 

The resources required for implementing Si can be deter- 
mined recursively. 



Ro(Si) 


= Ro{t e )^ 


-2 = 


4 


(34) 


Ri(Si) 


= Rl(te)-\ 


-3 = 


3 


(35) 


MSi) 


= R2(t e )-\ 


- 1 = 


2 


(36) 


Ro{Si+i) 


= 2R (Si) 


+ 2- 


f 2 1 - 1 


(37) 


Ri{Si+i) 


= 2Ri(Si) 


+ 4 




(38) 


R 2 {Si+i) 


= 2R 2 (Si) 


+ 1, 




(39) 



which one can solve to obtain 



Ro(Si) 


= 8x2' 


~ l -4 


(40) 


Ri(Si) 


= 7x2' 


- x -4 


(41) 


R 2 (Si) 


= 3x2' 




(42) 



The probability of success of Si+i using two independent out- 
puts of Si is (1 — g) 4 (l — p), which can be shown to imply 



polynomial total resource use. We now obtain new expres- 
sions for the Z k preparation resources (assuming that the last 
correcting series of Ziso's is deferred): 



7x2 



R (Z 2 ) = R (Si) + 1 
R^Z 2 ') = R^Si) 
R 2 (Z 2 ') = R 2 (Si) = 3 x 2 1 - 1 



i-i 



2 = 7 x 2 1 - 1 - 2 



1. 



(43) 
(44) 
(45) 



The success probability given the output of Si is (1 — q) 2 . 

To prepare the state needed for recovery after one of the 
qubits has been measured, follow the part of the protocol that 
does not involve the remaining qubit. There are two measure- 
ments that need to be made, and we note that the state can be 
used for completing th e pro tocol regardless of the outcome. 
As explained in Sect. VIA, the preparation can be decom- 
posed into making a copy of \t e ) and of an eigenstate of Z 
and then measuring XX, correcting the outcome with a Ziso 
if ncessary. Using the implementation of the XX measure- 
ment above and counting the two Y rotations needed to obtain 
the X operators from Z's in the couplings, we get 

Ro(r e ) = R (t e ) + 3 = 5 (46) 
Rx{r e ) = R 1 {t e ) + 2 = 2 (47) 
R 2 (r e ) = R 2 (t e ) + 2 = 3. (48) 

The success probability is (1 — q) 2 (l — p) 2 . 

B. Resources for operations 

To follow the resource usage through several levels of con- 
catenation, it is necessary to determine the maximum resource 
usage for each category of operations. First are the one qubit 
180° rotations, state preparations and measurements in the Z 
or Y basis. Of these, state preparation has the highest resource 
requirements for i?i and R 2 . i?o is highest for the 180° rota- 
tions. This gives the following estimates: 



Ro(0) 
Ri(0) 
R 2 (0) 



2 
3 
1. 



(49) 
(50) 
(51) 



The next category consists of the 90° Z or Y rotations. For 
the Y rotation, we conjugate a Z rotation by logical X 99 's. 
The resource analysis is complicated by the need for using 
recovery operations on partial failures. The expected num- 
ber of ZZ rotations that need to be retried after failure and 
successful recovery of the first qubit is 1/(1 — /(l — /)), as 
/(l — /) is the probability of failing and then successfully 
recovering. To get a better bound on the number of directly 
contributing operations, note that the ZZ couplings are imple- 
mented in such a way that if the source fails, the target is not 
touched. Thus teleportation steps that fail at the source of the 
coupling that precedes the measurements need not be counted 
except when estimating total resources. Thus, the expected 
number of contributing recovery operations can be bounded 
by 1/(1 - /(l - /)) - 1 = /(l - /)/(l - /(l - /)) for fail- 
ures at the first qubit and f(l - /)(! - f 2 /(l - /(l - /))) = 
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f(l - /) 2 /(l - f(l - /)) for failures at the second qubit. 
Thus 

i?o(l) = 2 + /(l-/)(2-/)/(l-/(l-/))Jfc(r e ) 

< 2 + 5/(l-/)(2-/)/(l-/(l-/)) (52) 
flx(l) = /(l-/)(2-/)/(l-/(l-/)))J? 1 (r e ) 

< 2/(l-/)(2-/)/(l-/(l-/)) (53) 
R 2 (l) = 1/(1-/(1-/)) 

+ /(l-/)(2-/)/(l-/(l-/)))i? 2 (r e ) 

< 1/(1-/(1-/)) 

+ 3/(l-/)(2-/)/(l-/(l-/)). (54) 

The final category has the 90° couplings. Except for at most 
4 Xgo rotations needed to get Y operators, it suffices to deter- 
mine the requirements for the logical ZZ operation. Let 2 l be 
the number of copies of \t e ) used in the prepared state. The 
calculation is similar to that for R x (l). Noting that the proba- 
bility of a teleportation failing at the target and the subsequent 
recovery succeeding is /(l — /) 2 , one can bound the expected 
number of recovery attempts by 4/(1 — /) 2 / (1 — /(l — /) 2 ). 
Some of the teleportations are attempted multiple times (just 
like the rotation is attempted multiple times in the previous 
case) and we need to account for the correction steps of the 
teleportation. 

Ro(2) = 8 + 2/(1- f(l-f))+Ro(Si) 

+ (4/(l-/) 2 /(l-/(l-/) 2 ))i? (r e ) 

< 8x2/(l-/(l-/)) + 2'- 1 

+ 20/(l-/) 2 /(l-/(l-/) 2 )+4 

(55) 

Ri(2) = 

+ (4/(l-/) 2 /(l-/(l-/) 2 ))i? 1 (r e ) 

< 2'- 1 + (8/(l-/) 2 /(l-/(l-/) 2 ))-2 

(56) 

R 2 (2) = 1/(1-/(1-/)) +^(5,) 

+ (4/(l-/) 2 /(l-/(l-/) 2 ))i? 2 (r e ) 

< 1/(1 -/(l-/)) + 3 X2 1 - 1 

+ 12/(1 -/) 2 /(l-/(l-/) 2 )-l. 

(57) 



and //(l - /(l - /)) = .308, .152. We evaluate the values 
for the coupling evolution at the second level. Using paren- 
thesized superscripts to denote the levels gives, for example, 

i?o (2) (2) = (7 x 2 2 + 20 x 0.156 + 9)i? (1) (0) 
+ (2 2 + 8 x 0.156 - 2)i? (1) (l) 
+ (3 x 2 2 + 12 x 0.156 + .5)i? (1) (2) 
< 656 (58) 

Similarly, 

i?i (2) (2) < 169 (59) 
i? 2 (2) (2) < 239 (60) 

Much of the inefficiency comes from the lack of optimization 
for implementing the logical ZZ rotation. 

The expected total resource usage including those needed 
for the failed state preparation attempts can be determined 
from the probabilities of successful preparation and requires 
recalculating the expressions for the various resources. For 
example, at the first level, p = pi = 0.44 and q = q\ = and 
at the second, p = P2 = -25 and q = qi = .135. Thus the ex- 
pected number of attempts required to make the state needed 
for error recovery is approximately 3.16 at the first level and 
2.38 at the second. 

The resource values imply that a controlled sign flip at the 
top level depends on less than 250 controlled sign flips imple- 
mented with pairs of three photon entanglements in iLOQC. 
This can be used to bound the effect of errors that cannot oth- 
erwise be controlled. The resource bounds improve if four or 
five photon states \t n ) can be reliably generated, thus giving 
better initial values of / and substantial gains in efficiency at 
the higher levels. 

VIII. COMPENSATING FOR OTHER ERRORS 

So far we have shown how to use a simple code with care- 
ful implementation of the basic operation to rapidly boost the 
probability of successful completion of gates. The next step is 
to consider the contribution of other errors and how to correct 
for them. Two important types of errors are phase errors and 
loss of particles. 



C. An explicit example 

Suppose the goal is to have a 5% failure probability per 
qubit, which is not far from the current best estimates for 
the communication threshold [|n|. If we use pairs of three- 
photon entanglements to generate the controlled sign flip at 
the first level in LOQC, the initial failure probability per 
qubit is / = 1/4. Using the expression for Fz in Eq. £3] 
gives / = 0.135 after one level and / = 0.038 after two. 
For simplicity, we choose I = 3 in both levels (hopefully a 
safe choice, though in principle this effects the probabilities 
a bit). Other useful values at the first and second levels are 
//(I - /) = 0.333, 0.156, 1/(1 - /(l - /)) = 1.23, 1.132 



A. Phase errors 

The occurrence of a phase error can be detected for the 
code used above by following the teleported XX measure- 
ment procedure in the absence of a failure. A — 1 eigenvalue 
indicates an error. At this point, it is necessary to return the 
state to the code space by applying an appropriate 180° rota- 
tion. Although it is not possible to correct for the error (we 
don't know which qubit is faulty), its detection can be used 
as information for correction in a higher level erasure code. 
An alternative is to use the generalization of the code to three 
qubits, which does have the capability of correcting for phase 
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errors. In fact, the fc-fold concatenation of the two qubit code 
with itself is actually a 2 fe qubit code that corrects for up to 
2 fc_1 — 1 phase errors. That property can be exploited by 
introducing an appropriate error detection/correction proce- 
dure periodically after the first few levels, without changing 
the overall concatenation scheme or the basic methods for im- 
plementing operations. 

B. Loss of particles and detector inefficiency 

The methods of LOQC include one that can, with some 
probability of failure, detect loss of a photon used to define 
a photonic qubit. The failure mode is again one involving a Z 
measurement, so the scheme can be used with the two qubit 
code. This will have an effect on the overall error behavior. 
We leave the calculations as an open problem. 

One observation not made in [[!]] is that if after any of 
the teleportation steps used for basic LOQC, we measure the 
modes not containing the teleported qubits, and the total num- 
ber of photons detected is not equal to the number initially 
prepared in the entanglement, then a photon was lost, possi- 
bly due to detector inefficiency. Such an event can be declared 
as a qubit loss. Doing so turns the problem of detector ineffi- 
ciency into one of having to handle detected qubit loss at some 
rate. It is a useful task to determine the relationship between 
detector inefficiency and the probability of detected loss. 

If total loss is detected, the two qubit code is insufficient 
for restoring the encoded information. Instead, it is necessary 
(perhaps after a few levels of two qubit encodings) to use era- 
sure codes. These are codes with the property that one (or 
a few) lost qubits can be restored without loss of informa- 
tion. The accuracy threshold for the error model where each 
operation satisfies that the target qubits are lost with indepen- 
dent probability s appears to be very good also, perhaps below 
95%. A brief explanation based on conservative calculations 
giving a value close to 99% is below. The erasure model is 
one of few for which it is possible to establish exact quantum 
communication quantities, such as the quantum channel ca- 
pacity [^8J. We therefore suggest that calculating the thresh- 
old for the erasure error-model is an excellent open problem 
to solve. 



C. An erasure code 

A useful erasure code encoding one qubit into four is de- 
fined by the stabilizer group generated by 

S E = {X^X^X^ 3 \Y^Y^Y^,Z^Z^Z^} (61) 

One can define encoded operators by 

Z (L) = Z (2) Z (3) ( 62 ) 

X W = x^X^. (63) 

We chose the erasure code and the logical operators for their 
small support (number of non-identity Pauli operators in the 



products). If it is desirable to encode two qubits into four, 
the erasure code with stabilizer generated by XXXX and 
ZZZZ can be used. 

It turns out that it is easier to analyze thresholds for erasure 
errors than for Z measurements. Implementations of opera- 
tions can be based on teleportation the same way we did be- 
fore, again relying on the ability to guarantee prepared states 
by retrying after failure. For the present discussion, we im- 
plement the teleportations required for an operation all in one 
step, and then follow up with error recovery if necessary. Of 
the operations required for quantum computing, state prepara- 
tion can be implemented error-free at all levels by retrying the 
process until success. To measure Z^ L \ first measure Z^ 2 \ 
and if this fails, restore the logical qubit. After successfully 
measuring , measure Z®\ and if that fails, and^ 4 '. 
The value of the Z^> measurement can be inferred by us- 
ing the parity constraint associated with the third generator of 
Se- The other logical Pauli operators can be measured sim- 
ilarly. The probability that the measurement fails is bounded 
by sr + 3s 2 , where r is the probability that the error recovery 
procedure fails. 

Suppose that qubit 1 is lost. The error recovery procedure 
requires measuring the first and third operators of Se, which, 
by reintroducing a fixed state for qubit 1, can be done with 
two full teleportation steps each. We ignore the fact that some 
failures in these steps can be recovered and estimate that each 
teleportation step fails with probability at most 3s (due to the 
coupling rotation and two measurements, not including the 
correction for now). The recovery step concludes with a cor- 
recting operation on qubit 1, so that we can estimate r < 13s 
(ignoring the possibility of retrying the recovery step if one of 
the last corrections fail). The probability of a failed measure- 
ment is therefore bounded by 16s 2 . 

Since we are unable to implement arbitrary rotations di- 
rectly, it is necessary to implement (say) a 45° rotation by 
teleportation. (The compensation step is now a 90° coupling 
rotation, which can be implemented either by teleportation or 
directly, exploiting the independence assumption on loss of 
qubits.) Since the most complex operation involves coupling 
two logical qubits with a 90° rotation, we finish by giving a 
rough estimate of this failure probability. The Hamiltonians 
to be evolved have weight four (two on each pair of encod- 
ing qubits), so it suffices to teleport four qubits after a suitable 
state preparation. If one of the teleportations fail, the recov- 
ery procedure is followed. We may need to compensate for 
a negative rotation induced by the teleportations' last correc- 
tion step, but that can be absorbed into the procedure. Thus 
the error probability for the first logical qubit is bounded by 
(8s)(13s) = 104s 2 , whence an accuracy threshold better than 
about 99%. 



D. General errors 

The techniques discussed so far compensate for all of the 
primary errors that are expected to occur in an LOQC sys- 
tem. Additional errors can be attributed to improper settings 
of beam splitters, undetected loss, stray photons, etc. The first 
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of these depends on accurate calibration of classical control 
parameters. Happily, this type of error affects the probabilities 
quadratically, so that it can be minimized well by engineering. 
To deal with other errors eventually requires more powerful 
quantum error correction, and the goal of any implementation 
is to minimize the need for these techniques. 

IX. CONCLUSION 

This work is a first attempt at reducing the overhead and 
need for efficient devices for implementing LOQC. It shows 
that even without much optimization of the operations or tight 
estimates of errors and resources, there are techniques that 
can be used to obtain useful quantum operations in LOQC 
with overheads within two orders of magnitude of those re- 
quired for other reliable implementations of quantum comput- 
ers. The advantages of LOQC include the ability to compen- 
sate readily for the primary errors in optics while preparing 
large numbers of the basic quantum systems, which are pho- 
tons in a superposition of two modes. It is necessary to further 
optimize the methods and to better analyze the error behavior, 
particularly the effects of detector inefficiencies, single pho- 



ton state preparation errors and timing or overlap problems 
for photons. A fruitful area of further investigation is to deter- 
mine whether larger codes and the method of encoding multi- 
ple qubits at once can be used to improve efficiency and error 
behavior. 

Our proposal consists of a multi-level system with various 
types of error-control gradually introduced and supported by 
high-output state preparation factories that exploit the ease 
with which many photons can be produced. Although in 
the very long term, solid state or molecular computing meth- 
ods are the preferred implementation for large scale quantum 
computing, LOQC is now a viable alternative to achieving 
the capability of non-trivial quantum information processing. 
One area where LOQC has a long term future is in commu- 
nication. Photon based systems are currently the only rea- 
sonable proposals for long distance quantum communication. 
Since the necessary accuracies for successfully exchanging 
entanglement over arbitrary distances are well below 99%, 
this may also be the first application of LOQC to quantum 
information processing to be experimentally implemented. 
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